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基于 高 亚洲 精细 再 分 析 数 据 模拟 普 若 六 日 
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摘 要 : 采用 高 亚洲 精细 再 分 析 数据 (High Asia Refined analysis, HAR ) L4 2% Coupled Snowpack and Ice Surface En- 
ergy and Mass balance Model (COSIMA ) 模 拟 了 青藏 高 原 最 大 冰原 一 一 普 知 岗 日 2012-04-23 一 2014-05-11 的 能 量 - 


物质 平衡 ， 


分 析 了 普 若 岗 日 冰原 消融 机 理 。 模 拟 的 2012-04-23 一 2013-05-23 与 2013-05-23 一 2014-05-11 的 物 


质 平衡 分 别 为 -0.18 m w.e.、-0.36 m w.e.。 在 整个 研究 时 段 内 ,冰原 西部 、 中 北部 等 海拔 低 于 5800 m 的 区 域 处 于 物质 
亏损 状态 ,最 大 达到 -2.69 m w.e.; 冰 原 北部 区 域 因 为 海拔 较 高 的 原因 ,处 于 物质 积累 状态 ,在 东北 部 区 域 最 高 值 达 


到 0.85 m w.e., 2013-05-23—2014-05-11 用 于 消融 的 能 量 相 比 2012-04-23 一 2013-05-23 减少 ,冰川 表面 消融 减 
速 ,但 普 若 岗 日 冰原 上 空 2013-05-23 一 2014-05-11 的 固态 降水 量 少 于 2012-04-23 一 2013-05-23 ,因此 冰川 在 


2013-05-23—2014-05-11 日 的 表面 


m 


在 气温 快速 上 升 的 背景 下 ,全 球 各 地 冰川 均 呈 
现 出 不 同 程度 的 消融 现象 一 ,位 于 青藏 高 原 内 部 的 
山岳 冰川 更 是 如 此 ”2 。 喜 马 拉 雅 山 附 近 的 冰川 物质 
损失 速率 在 2003 一 2008 年 间 达 到 了 0.59 m w.e.ar 
祁连山 地 区 的 冰川 物质 损失 速率 在 2003—2009 年 
间 达 到 了 0.77+0.35 m w.e.a , EDGE Jae JH 
地 区 的 冰川 D TE lares D RII E pa. H ok Jak cb 
于 比较 稳定 的 状态 。 已 有 一 些 学 者 通过 大 地 测量 
的 方法 对 普 知 岗 日 冰原 的 物质 平衡 变化 进行 了 定 
量 估计 。Neckel 等 "通过 差分 TanDEM-X 数 据 生成 
的 DEM 和 X 波 段 SRTM DEM ,估算 出 普 若 岗 日 冰原 
2000—2012 年 之 间 的 年 物质 平衡 值 为 -0.044+ 
0.015 m w.e.a!, Liu 等 串通 过 TanDEM-X 数 据 和 大 
地 测量 法 佑 算出 善 若 岗 日 冰原 在 2011 一 2012 4F 
2012—2013 年 2013—2014 年 .2014 一 2016 年 之 间 
的 物质 平衡 分 别 为 0.44+0.10 m w.e.a!,—0.13+0.03 
m w.e.a ! ,—0.3440.06 m w.e.a ' 以 及 -0.52+0.10 m w.e. 
a'。 大 地 测量 手段 受 限于 卫星 影像 数据 的 获取 , 物 


收 稿 日 期 : 2020-11-02; 修订 日 期 : 2021-04-28 


面 净 物 质 平衡 反而 减少 。 
关键 词 : MHRA; 高 亚洲 精细 再 分 析 数 据 ; COSIMA; 


能 量 平衡 ; 物质 平衡 


质 平衡 观测 时 间 分 辨 率 也 有 限 。 基 于 气象 参数 和 
辐射 值 的 能 量 - 物 质 平衡 模拟 方法 不 但 可 以 获取 精 
细 化 的 物质 平衡 估计 结果 ,还 可 以 从 物理 消融 机 制 
的 角度 剖析 冰川 变化 的 原因 ,预测 冰川 未 来 状态 ， 
更 深入 地 研究 冰川 对 气候 的 啊 应 关系 ”。 

冰川 表面 物质 平衡 主要 由 物质 积累 和 消融 构 
成 。 冰 川 消融 主要 受 控 于 冰川 表面 的 能 量 收 支 状 
况 ,冰川 积累 则 主要 来 源 于 外 部 降水 。 气 温 作为 直 
接 影响 冰川 物质 消融 的 重要 参数 , 它 在 冰川 物质 乎 
衡 模拟 中 扮演 着 重要 的 角色 。Hock ”提出 了 基于 
温度 变化 的 度 日 模型 ,用 于 模拟 冰川 物质 平衡 变 
化 。Huintjes 等 通过 引入 潜在 太阳 短波 辐射 对 度 
日 模型 进行 修正 ,但 改正 后 的 度 日 模型 仍然 只 涉及 
到 表面 的 物质 平衡 计算 ,并 没有 涉及 到 冰川 内 部 的 
消融 。Huintjes 等 "提出 Coupled Snowpack and Ice 
surface energy and Mass balance model (COSIMA ) 模 
型 。 该 模型 考虑 了 气压 、 温 度 、 风 速 等 冰川 表面 因 
素 以 及 表面 以 下 的 10 m 内 的 能 量 -物质 平衡 过 程 ， 
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消融 过 程 。 此 外 ,该 模型 的 时 间 分 辨 率 达 到 1h ,可 
以 精细 模拟 冰川 动态 变化 ,有 助 于 深度 揭示 冰川 能 
量 -物质 变化 过 程 。Huintjes 等 通过 该 模型 和 高 
亚洲 精细 再 分 析 数 据 (High Asia Refined analysis， 
HAR ) 估 算出 2001 一 2011 年 普 若 岗 日 冰原 的 消融 速 
率 为 -0.044 m w.e.a ', 与 Neckel 等 ”估计 的 结 
合 。 

普 知 岗 日 作为 青藏 高 原 中 腹地 区 最 大 的 冰原 ， 
是 气候 和 环境 变化 的 记录 者 ,同时 是 周围 湖泊 的 主 
要 供给 来 源 。Lei 等" 曾 利用 航拍 地 形 图 .SRTM 等 
数据 估算 普 若 岗 日 冰原 物质 平衡 变化 并 分 析 其 对 
令 臣 湖水 量 的 影响 ;Liu 等 利用 TanDEM-X 数 据 估 
算 了 普 若 岗 日 冰原 2012 一 2016 年 厚度 变化 。 但 这 
些 研 究 缺 乏 对 普 若 岗 日 冰原 近期 消融 机 理 的 分 
析 。 为 揭示 普 在 岗 日 近期 动态 ,进一步 分 析 其 对 气 
候 变化 的 响应 方式 ,本 文采 用 高 亚洲 精细 再 分 析 数 
据 (High Asia Refined analysis, HAR ) 4 F COSIMA 
模型 研究 了 2012—2014 E [8] 3&5; pg H vKU Be te AC 
CRUCE EE foL YES T PE Pd H VERUS 2012— 
2014 年 间 各 能 量 通 量 数据 。 本 研究 可 提供 高 时 间 
分 辨 率 的 普 春 疯 日 冰原 物质 平衡 变化 估计 结果 ,为 
分 析 普 铬 阅 日 冰原 物质 平衡 对 气候 变化 的 响应 方 


72°0'E 84°0'E 96°0'E 


88°36'E 88°48'E 


RARIS Ae TB AK Iz 28 ERU SE WU e BE Be Dt BESS 
考 数据 。 


1 数据 和 方法 


1.1 研究 区 概况 

普 若 岗 日 冰原 位 于 西藏 自治 区 双 湖 特别 区 的 
东北 部 ,覆盖 33?44'~34*04'N ,89°20'~89°50’E, Mil 
积 约 为 422.58 km"), PE) AE Kh ABO SH , 
最 高 海拔 为 6482 m"1;50 多 条 长 短 不 一 的 冰 舌 从 冰 
原 的 中 心 向 四 周 溢出 (图 1) , 冰 舌 最 低 处 海拔 约 
5350 m", E 1974 一 2010 年 间 冰 原平 均 雪 线 海拔 
高 度 约 为 5735 m5。 从 遥感 影像 上 看 , 普 若 网 日 冰 
原 表面 较为 洁净 , 冰 研 物 较 少 ,与 周围 地 物 容易 区 
分 开 来 。 普 车 岗 日 冰原 主要 受 印度 季风 和 中 纬度 
西风 带 控制 ,属于 极 大 陆 型 冰川 。 受 春季 西风 和 
夏季 印度 季风 的 影响 , 普 若 岗 日 峰 区 春季 和 夏季 降 
水 大 部 分 为 固态 降水 ,因此 善 若 岗 日 冰原 为 春 夏 混 
CURRUS 
1.2 数据 来 源 与 处 理 
1.2.1 高 亚洲 精细 再 分 析 数 据 高 亚洲 精细 再 分 析 
数据 (HAR ) 是 德国 柏林 工业 大 学 Dieter Scherer 教 
授 团 队 利用 数值 天 气 预 报 和 动力 降 尺 度 方法 获得 
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LE RAA RE RFE AL ERN R CE Ko COSIMA 模型 所 用 HAR 数据 区 域 (10 kmx10 km) 
图 1 38 pg Aa Ke A 


Fig. 1 Optical image of the Puruogangri ice cap 
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WAR BE”, EK MODIS 产品 以 及 气象 站 数 
据 进行 融合 ,结合 经 典 物理 参数 化 方案 ,采用 双向 钥 
套 模式 对 整个 亚洲 高 山区 进行 气象 参数 模拟 3。 数 
据 空间 分 辩 率 达 0.1°x0.1° ,时 间 分 辩 率 达 1h。 高 亚 
洲 地 区 数据 可 履 盖 2001 一 2014 年 间 ( 最 新 版 数据 已 
更 新 至 2018 年 )。 该 数据 精度 在 青藏 高 原 冰川 表面 
物质 平衡 研究 中 已 经 得 到 较 好 的 验证 55 。 本 研究 
下 载 了 普 若 岗 日 冰原 区 域 2012-04-23 一 2014-05- 
11 的 入 射 短波 辆 射 (W.m2) 2m 气 温 (K) .2 m UK 
(hPa) 10 m KUE (m* s) 、 降 水 (mm) 、 比 湿 、 云 量 以 
K SRTM DEM( 空 间 分 辨 率 :30 mx30m) 数 据 。 由 于 
二 者 分 辩 率 存在 差异 ,为 更 好 地 提取 普 若 岗 日 冰原 
范围 内 的 最 高 海拔 格 网 点 所 对 应 的 气象 数据 信息 ， 
方便 后 续 模 型 输入 数据 的 创建 ,本 研究 将 SRTM 
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DEM 重 采样 为 0.1°x0.1° 的 DEM ,并 基于 重 采样 后 的 
DEM 数据 提取 了 最 高 海拔 格 网 点 所 对 应 的 HAR 气 
象 数据 。 这 些 数据 随时 间 变 化 情况 见 图 2。 

1.2.2 数据 预 处 理 ”采用 Huintjes |" EHH pg 
日 冰原 2001 一 2011 年 间 物 质 平衡 时 所 使 用 的 梯度 
插值 方法 以 及 相应 气象 数据 的 梯度 因子 。 其 中 , 根 
据 HAR 数 据点 位 海拔 信息 与 对 应 点 位 的 比 湿 信息 ， 
对 所 有 点 位 进行 线性 拟 合 , 获 取 比 湿 梯度 因子 与 相 
关系 数 。 为 更 精确 地 体现 普 若 岗 日 冰原 空间 分 布 
上 的 消融 情况 ,本 研究 对 SRTM DEM( 空 间 分 辩 率 : 
30 mx30 m ) 进 行 重 采样 ,获取 空间 分 辨 率 为 450 mx 
450 m 的 新 SRTM DEM ,随后 以 最 高 海拔 格 网 点 为 
起 始点 ,结合 表 1 中 所 有 梯度 因子 新 SRTM DEM 以 
及 最 高 海拔 格 网 点 所 对 应 的 所 有 气象 数据 进行 冰 


T goo [ @) 短波 辐射 
E 300 
goo DEM TE ^ 0] 708 ME yee 
ae 
* EP a AMMAN Malen M 
oop WW. o, YAU al A ! " M PS 
TI 
5 500 Mo a p 
y 495 
490 E or op x ox i-jo ul p rc qo É a p cp dq 6. 45 1 5]. p ye uo cA su e i 
= 12r (a) 10m 风 速 
io i | 
3 1 fd ep MU A | MA AA a M j^ 
do 6 
a: 


气温 用 


降水 /mm 
E 


25 [ (e) 降水 量 


2012-04—08 2012-07-16 2012-10-23 2013-01-30 2013-05-09 2013-08-16 2013-11-23 2014-03-02 
时 间 / 年 -月 -日 


注 :降水 为 累积 值 外 ,其 余 均 为 日 均值 。 


图 2 普 若 岗 日 2012-04-23 一 2014-05-11 最 高 海拔 格 网 点 ( 选 自分 辩 率 为 0.1°x0.1° 的 DEM) 所 对 应 的 HAR 的 气象 数据 
Fig. 2 Meteorological data for the highest elevation grid point (selected from DEMs with a resolution of 0.1° x 0.1°) corresponding 


to the period 23 April 2012 to 11 May 2014 for Puruogangri 
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面 二 维 梯度 插值 ,最 终 获 取 冰 面 二 维 气 象 数 据 集 
(分 辨 率 .450 mx450 m)。 梯 度 因子 以 及 相关 系数 
如 表 1 所 示 。 

辐射 量 是 直接 影响 冰川 物质 消融 速率 的 重要 
参数 。 本 文 根 据 Kumar 等 在 1997 年 提出 的 辐射 
模型 ,结合 普 若 岗 日 地 理 位 置 海拔 、 纬 度 .坡度 以 
及 坡 向 等 信息 得 到 冰 面 晴空 人 射 短 波 辐射 (空间 分 
辩 率 为 450 mx450 m) ,根据 最 高 海拔 格 网 点 处 晴空 


表 1 HAR 气 象 数据 插值 所 采用 的 梯度 因子 
以 及 相关 系数 三 
Tab.1 Altitudinal gradients and correlation coefficients 
used to interpolate the grid HAR meteorological data 


(m); ww 为 10 m e BER) XE (m s) 6 

另外 , 湿 空 气 中 水 汽 重 量 与 于 空气 重量 的 混合 
比 会 影响 水 汽 转换 过 程 ,进而 影响 潜 热 通 量 的 计 
算 。 我 们 利用 混合 比 与 比 湿 的 关系 ( 式 42) 得 到 混 
合 比 。 水 汽 压 可 以 影响 下 行 长 波 辐射 ,进而 影响 净 
长 波 辐射 的 计算 。 得 到 混合 比 后 ,我们 可 以 通过 式 
S" 计算 水 汽 压 。 


662q 


~ 622- 0.622q (4) 
_ M, | 0.662e 
Log gg ers (5) 


AP: q 为 比 湿 ; W IRAE; M, 为 湿 空 气质 量 
(g); M, 为 干 空气 质量 (g); e 为 水 汽 压 (hPa); p 为 


梯度 因子 R 
气温 /(K:m') —0.00830 0.98 
UR/(hPa m") -0.06700 1.00 
KEJIK/(%-m™) 0.05300 0.67 
VERE (Io m) -0.01050 0.96 
10 m 风速 /ms m^) 0.00170 0.13 
zi 0.00001 0.09 


入 射 短波 辐射 值 与 HAR 入 射 短波 辐射 值 计算 转换 


因子 ,从 而 获得 冰 面 HAR( 空 间 分 辨 率 为 450 mx 
450 m) 人 射 短 波 辐射 值 。 
i, = Woss (1) 
"o SW na mn HAn 


式 中 : ri,, WERIT; SWW, ,为 对 应 点 位 的 晴空 辐 
射 值 ; SW, asas 为 对 应 点 位 的 HAR 辐射 值 。 

固态 降水 是 影响 冰川 物质 积累 的 关键 参数 。 
Miller 等 “在 2007 年 提出 了 双 曲 正切 函数 。 该 函数 
可 以 实现 +1 %C~+5 % 之 间 的 固态 和 液态 降水 分 离 。 
基于 此 函数 和 HAR 降水 数据 获得 了 普 知 岗 日 冰原 
的 固态 降水 量 。 公 式 如 下 : 

snowg — (2) 

prep (0.5 (-Ctanh((T2 - 273.15) - 2.5) x 2.5) - 1) 
式 中 : snowg 为 固态 降水 量 (m); prep 为 总 降水 量 
(m); T2 为 2m 的 温度 (K)。 

风速 湿度 等 参数 影响 着 冰川 表面 的 清流 交 
换 , 对 冰川 表面 能 量 交 换 过 程 有 一 定 的 影响 。HAR 
数据 集 只 提供 了 10 m 高 度 的 风速 。 本 人 研究 采用 Al- 
len 等 “在 1998 年 提出 的 风速 转换 模型 进行 不 同 高 
度 的 风速 转换 : 

u, = (u,,4.87)/[1og(67.8z — 5.42)] (3) 

IP: w, 292 m i REB UE (m «80 5 2 Zr ESL e RE 


ES 


气压 (hPa)。 
初始 雪 深 影响 模型 中 地 下 前 面 层 的 划分 。 对 
于 初始 雪 深 的 确定 ,采用 Huintjes 等 所 模拟 普 若 岗 
日 冰原 表面 能 量 -物质 平衡 时 采用 的 方法 ,即将 物 
质 平 衡 线 以 下 区 域 的 初始 雪 深 设 为 0 em, 物质 平衡 
线 以 上 至 海拔 最 高 点 区 域 的 初始 雪 深 通过 线性 插 
值 估计 。 
1.3 研究 方法 
通过 COSIMA 模型 和 气象 数据 模拟 冰川 的 消融 
速率 ,再 结合 固 液 态 降 水 模拟 冰川 物质 积累 ,从 而 
得 到 冰川 物质 平衡 ,分 析 冰 川 质 量变 化 对 气候 变化 
的 响应 。 
COSIMA 模型 对 冰川 的 能 量 和 物质 平衡 进行 分 
层 分 析 。 其 能 量 平衡 模型 如 下 : 
PFSBW. + LW GOL + Qa + Oe (6) 
式 中 : F OAKS BRE AE Ht 5 SW, 为 净 短波 辐射 量 ; 
LW 为 净 长 波 辐射 量 ; Qo 为 潜 热 通 量 ; 0, 为 感 热 
通 量 ; Q, 为 地 下 导热 量 。 
感 热 潜 热 的 计算 如 下 : 
Qon EPa et um. T) (7) 
Qu pa bas! Cu! ur (qu q.) (8) 
RP: p 为 空气 密度 ; 6c, 为 空气 等 压 比热容 
[1004.67 J* (kg-*C) ]; Ta . T, 分别 为 2 m 处 温度 以 
及 冰川 表面 温度 ; Lus NZS AAPA (Q.514Xx10* J kg / 
THEE IA (Q.849x105J-kg 05; qq. 2230272 m Ab 
比 湿 以 及 冰川 表面 比 湿 ; C, . Co 分 别 为 潜 热 块 体 
交换 系数 以 及 感 热 块 体 交 换 系数 (由 于 两 者 差别 小 
于 5%, 所 以 一 般 假设 两 者 相等 )。 块 体 交 换 系 数 可 
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通过 Braithwaite 在 1995 年 提出 的 公式 计算 : 


— k 
C= [In(z/z;) 9) 


式 中 :大 为 van Karman 系数 (0.41); z HWM se os es 
度 ; z 为 粗糙 长 度 。 

地 下 导热 量 主要 由 透射 短波 辐射 与 地 下 传 热 
构成 。 地 下 传 热 主要 根据 表面 与 次 表面 温度 差异 以 


及 传 热 系数 计 算 。 传 热 系数 的 确定 采用 Anderson'”” 
在 1976 年 所 提出 的 方法 : 

A =0.021 + 2.5(1000"')’ (10) 
式 中 : 为 传 热 系数 ; p 为 雪 层 密度 。 


透射 短波 辐射 根据 Bintanja 等 在 1995 年 提出 

的 模型 估算 : 
S(ZJ=S 肥 (1-2)e” (11) 

WP: S, 为 短波 辐射 在 深度 Z 处 的 剩余 部 分 ; SW... 
为 冰川 表面 吸收 的 净 短 波 辐射 ; & 为 净 短 波 辆 射 在 
顶层 模型 中 被 吸收 的 比例 ( 冰 / 雪 分 别 为 :0.8/0.9); B 
为 消光 系数 ( 冰 / 雪 分 别 为 :2.5/17.1)。 

反照 率 决 定 着 冰川 表面 吸收 太阳 短波 的 辐 映 
量 。 模 型 中 反照 率 的 确定 采用 的 是 Oerlemans 等 
在 1998 年 提出 的 方法 : 


Qo =Q + (Ca 一 a, Jexp[ 2 (12) 


Gau FOO, exp c3 (13) 


式 中 :i 为 模型 开始 运行 到 第 一 次 降雪 的 时 间 (d); 
a! 为 第 i 天 的 表面 反照 率 ; as 为 第 i 天 的 积 雪 反 
BR; am 为 粒 雪 反 照 率 (0.55); ou... 为 新 雪 反 照 
3€(0.9) ; s 为 距离 上 次 降雪 时 间 ; 三 为 新 雪 变 为 粒 
雪 的 时 间 尺 度 (6 d); al 为 冰 面 反照 率 (0.3); d 为 
雪 深 ; d* 为 雪 深 特征 尺度 (8 em) 。 

采用 Huintjes 等 模拟 普 若 岗 日 冰原 能 量 - 物 
质 平衡 时 所 采用 的 分 层 方法 ,以 0.1 m 为 深度 将 地 
下 10m 齐 分 为 100 个 剖面 层 , 并 分 别 对 其 密度 il 
度 等 参数 进行 分 析 和 计算 。 采 用 Milg 等 ”在 2012 
年 提出 的 参数 化 方案 确定 初始 的 新 雪 、 粒 雪 、 冰 层 
的 反照 率 以 及 粗糙 长 度 。 初 始 底层 温度 则 沿用 
Huintjes 等 ”模拟 普 知 岗 日 冰原 能 量 -物质 平 衡 时 
所 采用 的 数值 。 在 整个 模拟 过 程 中 ,COSIMA 模型 
所 采用 的 参数 以 及 分 层 方法 均 已 经 过 Huintjes Sm 
rand, ETE Ar Pd H VKSUBUSHB E BU AC. 


ChinaXiv 合 作 期 刊 


4 期 李 建 江 等 ;基于 高 亚洲 精细 再 分 析 数 据 模拟 普 若 岗 日 冰原 2012 一 2014 年 表面 能 量 -物质 平衡 923 


2 结果 与 分 析 


2.1 表面 能 量 平衡 

获取 的 2012 一 2014 年间 普 阁 岗 日 冰原 表面 主 
要 能 量 通 量 如 图 3 所 示 。2012 一 2013 年 间 净 短波 
辐射 量 为 63.81 W.m”, 感 热 通 量 为 35.59 W ^m^, 1$ 
长 波 辐射 量 为 -50.04 W -m? , 潜 热 通 量 为 -32.26 W - 
m”, 地 下 传 热 为 -10.04 W*m?, 2013—2014 4F [8] 4 
短波 辐射 量 为 67.09 W «m^? , 感 热 通 量 为 36.54 Wem”, 
净 长 波 辐射 量 为 -53.82 W - m", RP Tit 7g -33.34 
W.nm2, 地 下 传 热 为 -10.31 W.m?。 将 各 能 量 通 量 结 
果 ( 图 3) 与 Huintjes 55 ES TIL BS 35 25; pg H 2001— 
2011 年 表面 各 能 量 通 量 结果 对 比 , 发 现在 2001 一 
2014 年 间 净 短波 辐射 量 \ 净 长 波 辐射 量 、 感 热 通 量 、 
潜 热 通 量 等 能 量变 化 趋势 一 臻 。 在 整个 研究 期 间 ， 
净 短 波 辐射 通 量 值 呈现 一 定 的 规律 ,分 别 在 在 冬 / 夏 
季 达 到 最 低 /高 峰 。 一 方面 可 能 是 由 于 地 球 自转 、 公 
转 使 得 普 若 岗 日 冰川 表面 在 冬 /夏季 接收 太阳 光照 
时 间 有 所 差异 ; 男 一 方面 可 能 是 因为 冬 /夏季 发 生 固 
态 降 水 事件 的 概率 差异 ,导致 冬 / 夏 季 冰 川 表面 反照 
率 差异 。 净 长 流 辐射 也 具有 一 定 的 规律 ,在 夏季 达 
到 最 低 峰 , 这 可 能 是 由 于 云 量 以 及 比 湿 在 夏季 增加 
(图 2b、 图 2e), 导 致 下 行 长 波 辐射 增 大 。 潜 热 通 量 
在 干燥 的 条 件 下 一 般 为 负 值 ; 在 夏季 时 ,风速 较 低 
(图 24), 不 利于 冰川 表面 与 空气 界面 之 间 的 注 流 交 
换 , 此 时 , 潜 热 通 量 为 正 值 ,冰川 表面 发 生 凝 结 或 凝 
华 现象 。 空 气 界面 与 冰川 表面 存在 气温 差 ,导致 空 
气 界 面向 冰川 表面 进行 能 量 传递 ,从 而 产生 感 热 通 
量 , 感 热 通 量 一 般 是 为 正 值 旦 大 小 与 潜 热 通 量 值 基 
本 持平 。 

表 2 列 出 了 2012-04-23 一 2014-05-11 各 能 量 
通 量 占 比 ,注意 能 量 总 和 为 各 通 量 绝对 值 之 和 。 相 
比 于 2012-04-23 一 2013-05-23 ,2013-05-23—2014- 
05-11 净 长 波 辐射 增加 3.78 W «m? , 净 短 波 辐射 增加 
3.28 Wm :, 感 热 通 量 增加 0.95 Wm”, HRON EIA 
加 1.08 W.m 2?, 地 下 传 热 增加 0.27 Wem, BEAK 
看 ,2013-05-23 一 2014-05-11 可 用 于 物质 消融 的 能 
量 减 少 。2012-04-23 一 2014-05-11 冰川 能 量 收 入 
主要 来 自 于 太阳 短波 辐射 (32.8% ) 与 感 热 通 量 
(18.1%), 而 能 量 支 出 主要 是 净 长 波 辐射 (26.1%) 与 
潜 热 通 量 (16.5%)。 净 长 波 辐射 增加 可 能 是 由 于 云 
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图 3 2012-04-23—2014-05-11 普 若 岗 日 15 d 各 能 量 通 量 均 值 
Fig. 3 Different components of the 15 days averaged energy fluxes at Puruogangri ice cap between April 23, 2012 and May 11, 2014 


表 2 2012-04-23 一 2014-05-11 普 若 岗 日 各 能 量 通 量 年 均值 能 量变 化 及 占 比 
Tab.2 Annual variation and contributions of different energy components at the puruogangri ice cap 
between April 23, 2012 and May 11, 2014 


量 的 适当 减少 (使 得 下 行 长 波 辐 射 减少 ) , 净 短 波 辆 
射 增加 可 能 是 由 于 固态 降水 减少 (使 得 冰 面 反照 率 


能 量 绝对 值 总 和 = MEE = 净 长 波 辐射 ” 感 热 通 量 ， 潜 热 通 量 ”地 下 传 热 ”消融 能 量 
2012-04-23 一 能 量 /(W.m”) 195.52 63.81 -50.04 35.59 -32.26 -10.04 -3.78 
2013-05-23 占 比 /% = 32.6 25.6 18.2 16.5 5.1 2.0 
2013-05-23— —— &ER/(- m?) 203.07 67.09 -53.82 36.54 -33.34 -10.31 -1.97 
2014-05-11 占 比 /% = 33.0 26.5 18.0 16.4 5.0 14 


#3 2012-04-23—2014-05-11 普 若 岗 日 物质 平衡 
分 量 模拟 值 


Tab.3 Different components of the simulated mass 


降低 )。 而 感 热 通 量 、 潜 热 通 量 以 及 地 下 传 热 变化 


量 较 小 的 原因 则 可 能 是 温度 、 气 压 、 比 湿 以 及 风速 
的 变化 较 小 。 


2.2 物质 平衡 变化 DEFER O23 20M 05 IL 
表 3 和 表 4 分 别 列 出 了 2012-2014 年 间 普 若 岗 PUTENT T ET 

日 冰原 物质 平衡 各 分 量 与 驱动 COSIMA 模型 所 需 各 冰 下 消融 /(m w.e.) -0.09 -0.13 

要 素 的 年 均值 /累积 值 。 相 比 2012 一 2013 年 ， 升华 /(m w.e.) -0.01 -0.02 

2013 一 2014 年 表面 消融 量 减少 , 冰 下 消融 量 增加 ， 融 水 再 冻结 /(m w.e.) 40.02 +0.02 

物质 升华 量 增加 , 融 水 再 冻 量 略微 减少 ,固态 降水 固态 降水 /Com w.e.) +0.57 +0.33 

减少 。 冰 川 表面 消融 减少 是 由 于 表面 消融 能 量 减 物质 平衡 /(m w.e.) -0.18 70.36 


balance of puruogangri ice cap between April 23, 2012 


and May 11, 2014 
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#4 2012-04-23—2014-05-11HAR 气象 数据 值 
Tab.4 HAR meteorological data between April 23, 
2012 and May 11, 2014 


参数 2012-04-23 一 2013-05-23 一 
2013-05-23 2014-05-11 
问 下 短波 辐射 /(W.m”) 256.5 245.6 
气温 /k 261.7 260.9 
气压 /hPa 500.9 501.3 
降水 量 /mm 752.8 639.9 
比 湿 /(g* kg) 2.9 2.8 
风速 /(ms') 6.9 6.7 
云 量 /% 30.8 29.0 


少 。 净 短波 辐射 在 向 下 透射 时 会 导致 表层 以 下 每 
个 剖面 层 出 现 消融 。 冰 下 消融 量 增加 是 由 透射 短 
波 辐射 增加 导致 。 物 质 升华 量 增加 由 潜 热 通 量 增 
加 导致 。 融 水 再 冻 量 略微 减少 是 因为 温度 的 略微 
下 降 , 造 成 冰川 表面 融 水 减少 。 固 态 降 水 减少 ,一 
方面 可 能 是 由 于 2013 一 2014 年 风速 下 降 ( 表 4)。 青 
藏 高 原 受 印度 季风 .中 纬度 西风 控制 ,而 印度 季风 
和 中 纬度 西风 将 南 印度 洋 .孟加拉湾 上 空 水 汽 输 送 

青藏 高 原 ””。 风 速 下 降 意味 着 输送 的 水 汽 量 减 
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17 号 冰川 整体 海拔 偏 低 ,3、4、11 号 冰川 的 冰 舌 海拔 
偏 低 , 但 上 部 海拔 较 高 。 海 拔 较 低 的 区 域 ,固态 降 
水 一 般 相 对 较 少 ,物质 累积 较 少 ; 冰 面 反照 率 偏 低 ， 
吸收 的 短波 辐射 量 偏 多 ,同时 空气 温度 相对 偏 高 ， 
感 热 通 量 偏 高 ,因此 用 于 消融 的 净 能 量 相对 较 高 ， 
冰雪 消融 更 严重 。 反 之 ,海拔 较 高 的 区 域 ,固态 降 
水 一 般 相对 更 多 ,物质 累积 更 多 ; 冰 面 反照 率 偏 高 ， 
吸收 的 短波 辐射 量 偏 少 ,同时 空气 温度 相对 偏 低 ， 
感 热 通 量 偏 低 , 用 于 消融 的 净 能 量 相 对 较 低 , 冰 雪 
消融 更 少 。 
2.3 不 确定 性 分 析 

COSIMA 模型 中 含有 未 知 的 参数 。 由 于 在 普 若 
岗 日 冰原 缺乏 实地 观测 值 ,无 法 对 模型 参数 进行 具 
体 评估 。 文 中 这 些 参数 的 初始 值 来 自 Malg 2E 
结 的 经 验 值 或 者 扎 当 冰川 的 实测 数据 "5。Huintjes 
等 所 采用 同样 的 参数 模拟 普 若 岗 日 冰原 2001 一 
2011 年 能 量 -物质 平衡 ,取得 较 好 的 效果 。Molg*”" 
等 在 2012 年 和 2014 年 对 模拟 类 似 于 扎 当 冰川 结构 
的 单 点 能 量 -物质 平衡 模型 进行 了 敏感 性 评估 ,其 
通过 改变 模型 的 29 个 参数 ,对 位 于 扎 当 冰川 的 自动 


少 。 另 一 方面 可 能 由 于 气温 下 降 ( 表 4) ,青藏 高 原 
内 部 湖泊 蔡 发 量 减 少 ,空气 中 水 汽 含量 相应 减少 ， 
此 外 , 表 4 中 比 湿 数据 也 可 看 出 ,空气 中 水 汽 含量 的 
减少 。 在 2012 一 2014 年 ,物质 积累 主要 来 源 于 固态 
降水 (0.90 m w.e.) ,其 次 是 融 水 再 冻结 (0.04 m w. 
e.)。 物 质 亏 损 主 要 由 表面 消融 (-1.23 m wee.) UF 
消融 (-0.22 m w.e.) 以 及 物质 升华 (-0.03 m w.e.) 
引起 。 

由 图 4 可 以 看 出 , 普 若 岗 日 冰原 西部 区 域 大 多 
处 于 物质 亏损 状态 ,最 大 损失 量 达 到 2.69 m w.e.; 普 
若 岗 日 冰原 东部 区 域 大 多 处 于 物质 积累 状态 ,最 大 
达到 0.85 m w.e.。 图 4 中 1~17 号 冰川 为 本 区 域 面 积 
KF 10 km 的 冰川 。 上 有 具体 来 说 ,1、3、4、11、14、15、17 
号 冰川 的 冰 舌 部 分 物质 亏损 较 严 重 ;2、5、6、7、8.9、 
16 号 冰川 上 部 出 现 物质 积累 。 整 体 的 物质 变化 空 
间 分 布 趋势 与 以 往 能 量 - 物 质 平衡 结果 趋势 较为 一 
致 引 。 表 5 展示 了 1~17 号 的 冰川 在 2012 一 2014 年 
的 物质 平衡 值 。1、3、4、11、14、15、17 号 冰川 物质 亏 
员 量 高 于 整个 冰原 平均 值 (-0.54 m w.e.a ') ,而 2、5、 
6.7.8.9、10、12、13、16 号 冰川 物质 亏损 量 低 于 整个 冰 
原平 均值 。1 号 冰川 物质 亏损 量 最 大 (-0.87 m w.e.),8 
号 冰川 物质 亏损 量 最 低 (-0.02 m w.e.) 4 1.14415, 


气象 站 所 观测 的 数据 进行 蒙特 卡 罗 模 拟 ,实验 证 明 
忽略 消融 期 积 雪 厚度 导致 的 误差 最 大 。 此 外 ,Milg 
等 3 通过 停 用 模型 中 部 分 物理 过 程 (短波 辐射 的 透 
射 , 满 流通 量 的 稳定 度 校 正 、 雪 的 压 实 和 沉淀 MV e 
层 的 融化 以 及 粗糙 长 度 的 变化 等 ) 分 析 物 质 平衡 模 
拟 结果 变化 ,发现 滑 流通 量 的 稳定 度 校 正 . 雪 的 压 
实 和 沉淀 . 亚 表 层 的 融化 对 物质 平衡 模拟 结果 影响 
较 大 。Huintjes EUF] COSIMA 模型 模拟 了 扎 当 
冰川 表面 能 量 -物质 平衡 结果 并 与 实测 结果 进行 对 
比 ,得 出 物质 平衡 结果 不 确定 性 为 600 kee m+ a's 
由 于 缺乏 实测 数据 ,本 文 无 法 定量 估计 模拟 结果 的 
不 确定 性 。 但 本 文 模拟 的 普 若 岗 日 冰原 2012 一 
2014 年 物质 平衡 定量 结果 和 空间 分 布 趋势 (图 4) 与 
Liu 等 ”通过 大 地 测量 方法 估计 的 结果 十 分 吻合 ， 
结果 可 靠 性 较 高 。 本 文 模拟 的 2012 一 2013 年 、 
2013—2014 年 累积 物质 平衡 分 别 为 -0.18 m w.e. 
和 -0.36 m w.e., Liu 等 中 估计 的 2012—2013 年 、 
2013—2014 年间 累积 物质 平衡 分 别 为 -0.13+0.03 m 
w.e. 和 -0.34+0.06 m w.e.。 通 过 COSIMA 模型 模拟 的 
普 若 岗 日 冰原 物质 平衡 值 与 Liu 等 加 的 定量 估计 结 
果 十 分 吻合 , 且 消 融 与 积累 的 空间 分 布 趋势 大 体 一 
致 ,因此 结果 可 靠 性 较 高 。 
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图 4 基于 COSIMA 模 型 模拟 2012-04-23 一 2014-05-11 普 若 岗 日 冰原 物质 
平衡 空间 分 布 结果 (空间 分 辨 率 :450 mx450 m) 


Fig.4 Distribution of the mass balance of the Puruogangri ice cap between April 23, 2012 and May 11, 


2014 simulated by COSIMA model (spatial resolution: 450 mx450 m) 


#5 普 若 岗 日 冰川 2012-04-23 一 2014-0S-11 的 物质 平衡 值 


Tab. 5 Cumulative mass balance of 17 glaciers within the Puruogangri ice cap between April 23, 2012 and May 11, 2014 


冰川 编号 01 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 
ane -0.87 -0.37 -0.59 -058 -0.30 -0.43 -0.18 -0.02 -0.27 -0.44 -0.76 -0.45 -0.48 -0.71 -0.67 -0.48 -0.73 
m w.e. 


此 外 ,Molg 等 ”将 HAR 数 据 中 的 降水 量 与 扎 当 
冰川 实测 的 降水 量 进行 了 对 比 ,发 现 HAR 降水 量 存 
在 过 高 估计 , 且 求 得 比例 因子 为 0.56。Huintjes 等 
采用 3 个 不 同 的 HAR 降水 量 比例 因子 (0.31.0.56、 
0.81) 分 别 进行 了 整个 普 知 岗 日 冰原 的 能 量 -物质 平 
衡 模拟 ,物质 平衡 结果 分 别 为 -2298 kg.m a '、-44 
kg:m?*a,*162 kgm2.a ,与 Neckel 等 ”获取 的 大 
地 测量 结果 (-44+15 kg.m?.a0) 相 比 可 知 ,HAR 降 


水 比例 因子 的 选取 对 模拟 结果 有 很 大 的 影响 , 当 比 
例 因 子 选 取 0.56 时 ,HAR 降水 数据 在 普 若 岗 日 冰原 
所 带 来 的 不 确定 性 较 低 。 


3 结论 


(1) 模拟 的 普 知 岗 日 冰原 在 2012- 04- 23 一 
2013-05-23 与 2013-05-23 一 2014-05-11 物质 平衡 
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他 学 者 
结果 (-0.13+0.03 m w.e.、-0.34+ 
0.06 m w.e.) 基 本 吻合 , 且 消 融和 积累 空间 分 布 趋势 
也 与 大 地 测量 结果 基本 吻合 ,结果 可 靠 性 较 高 。 
(2) 普 若 岗 日 冰原 主要 能 量 收入 项 为 净 短波 辐 
射 (32. OO aR 出 项 为 
净 长 波 辐射 (26.1% ) RITE BERE (16.596) HF 
2012-04-23 一 2013-05-23 ,2013-05-23—2014-05- 
11 的 净 短 波 辐射 变化 量 大 于 感 热 通 量 , 净 长 波 辐 身 
变化 量 大 于 潜 热 通 量 , 最 终 用 于 冰川 消融 的 能 量 
(3) 冰川 物质 累积 主要 来 源 于 固态 降水 ,其 次 
是 融 水 再 冻结 和 水 汽 凝 华 ,而 物质 亏损 主要 由 表面 
消融 和 冰 下 消融 构成 ,其 次 是 升华 作用 。 相 较 于 
2012-04-23—2013-05-23 ,2013-05-23—2014-05- 
11 的 冰 下 消融 量 增 加 ,表面 消融 量 减少 ,整体 消融 


分 别 为 -0.18 m w.e.,-0.36 m w.e., 与 早 前 其 
发 表 的 大 地 测量 


量 减少 ,但 由 于 固态 降水 量 的 减少 ,物质 净 损 失 量 
增加 。 
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Simulating surface energy and mass balance of the Puruogangri ice 
cap during 2012-2014 based on High Asia Refined analysis data 


LI Jianjiang, LI Jia, WU Lixin, WANG Yingzheng", GUO Lei' 
(1. School of Geosciences and Info-physics, Central South University, Changsha 410083, Hunan, China; 
2. College of Earth and Environmental Sciences, Lanzhou University, Lanzhou 730000, Gansu, China) 


Abstract: Glacier mass balance is an important parameter for describing how glaciers respond to climate change. 
The Puruogangri ice cap is the largest ice cap on the Tibetan Plateau, and the retrieval of its mass balance is of 
great significance. In this study, the surface energy balance and mass balance of the Puruogangri ice cap between 
April 23, 2012 and May 11, 2014 were simulated from the High Asia Refined analysis data via the coupled 
snowpack and ice surface energy and mass balance model(COSIMA). Based on simulated surface energy balance 
and mass balance, we revealed the ablation mechanism of the Puruogangri ice cap and the reason why ablation 
accelerated between April 23, 2012 and May 11, 2014. From April 23, 2012 to May 23, 2013 and May 23, 2013 
to May 11, 2014, the simulated mass balance values of the Puruogangri ice cap were —0.18 m w. e. and 0. 36 m w. 
e. , respectively. These results were consistent with the previously reported result (—0.13 + 0.03 m w. e. and -0. 34 + 
0.06 m w.e., respectively) that was derived by differencing the digital elevation models obtained at different 
times. The simulated mass balance indicated that the low areas, such as the western parts of the ice cap, 
experienced mass loss, while the high areas, such as the northern part of the ice cap, experienced mass gain. The 
maximum mass loss and gain during the observation period were 2.69 m w. e. and 0.85 m w. e. , respectively. The 
energy used for ablation between May 23, 2013 and May 11, 2014 was less than that used between April 23, 2012 
and May 23, 2013. Therefore, the surface ablation was weakened between May 23, 2013 and May 11, 2014. 
However, the solid precipitation between May 23, 2013 and May 11, 2014 was also less than that between April 
23, 2012 and May 23, 2013; therefore, the net surface mass balance of the ice cap decreased between May 23, 
2013 and May 11, 2014. 

Keywords: Puruogangri ice cap; High Asia Refined analysis; coupled snowpack and ice surface energy and mass 


balance model(COSIMA); surface energy balance; mass balance 


